Photon-graviton pair conversion 
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We consider the conversion of gravitons and photons as a four-wave mixing process. A nonlinear 
coupled systems of equations involving two gravitons and two photons is obtained, and the energy 
exchange between the different degrees of freedom is found. The scattering amplitudes are obtained, 
from which a crossection for incoherent processes can be found. An analytical example is given, and 
applications to the early Universe are discussed. 
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The interaction between electromagnetic waves 
(EMWs) and gravitational waves (GWs) have been 
considered by many authors, using approximate meth- 
ods 0, or by deriving exact solutions 
0. The later solutions are usually found on a flat 
vacuum background, in which case the simplest non- 
linear interaction processes involving three waves are 
excluded due to energy-momentum conservation. In 
the present letter we will consider four-wave interaction 
on a Minkowski background, the lowest order process 
allowing for energy exchange between photons and 
gravitons. Using the classical Einstein-Maxwell system, 
we adopt a calculation scheme inspired by quantum 
mechanics GH El E3 Assuming that two GWs 
interact with two EMWs, the only combination leading 
to energy exchange, the coupled system of equations is 
solved perturbatively, taking into account terms up to 
third order in the amplitudes. The four wave interaction 
coefficients (giving the quantum mechanical transition 
rates) are calculated in the center of mass system for 
arbitrary polarizations and scattering angle. Thus, the 
local energy transfer rate between GWs and EMWs is 
found, i.e. the probability of photoproduction by two 
gravitons or that of the inverse process. 

Although the time-scale for coherent GW-EMW pair 
conversion is shorter than the incoherent one, the large 
energy densities in the early universe makes it interest- 
ing to consider incoherent phenomena involving multi- 
ple conversion processes. Using the analogy between 
the interaction coefficients and the scattering matrix el- 
ements, the differential cross-section for graviton-photon 
pair conversion is calculated. In particular, it is found 
that a thermalization between photons and gravitons due 
to pair conversion requires a temperature T — > T p , where 
T p is the Planck temperature. 

In this letter, we use metric signature ( — h++), tetrad 
indices a, b, . . . = 0, 1, 2, 3, a, j3, . . . = 1, 2, 3, and coordi- 
nate indices fi, v, . . . = 0,1, 2, 3. The velocity of light in 
vacuum is c = 1. 

We use the formalism developed in Refs. 0,000]; 
introducing an observer four- velocity, V a , and decom- 
posing the electromagnetic field into an electric and 
a magnetic part. Introducing an orthonormal frame 



{e a - e£0 M }, where e = V = V a e a , (i.e. V a = 8$), 
the Maxwell field equations can be written in terms of 
the electric and magnetic fields with the inclusion of ef- 
fective charges and currents originating from the gravi- 
tational field (see e.g. 000). 

From the field equations one may derive the following 
inhomogeneous wave equations 



UE a = -e j E 
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where □ = eoeo — V ■ V and where the three- vector nota- 
tion E = (E a ) = (E 1 .^ 2 .^ 3 ) etc. and V = (ei,e 2 ,e 3 ) 
from 000 is used. The C§ c are the commutation 
functions, i.e. [e„,e&] = C^ b e c ). 

We consider plane waves with weakly space and time 
dependent amplitudes, i.e. E = E(x^) exp [ifc^a^] + c.c. 
and \d fl E(x l/ )\ -C |fc M ||i?(a; ly )|, where c.c. denotes the 
complex conjugate. In a flat vacuum background, energy- 
momentum conservation excludes three wave processes 
leading to energy exchange between photons and gravi- 
tons. Furthermore, the only four-wave combination con- 
verting energy between the gravitational and electromag- 
netic degrees of freedom involves two EMWs and two 
GWs. We denote the gravitational and electromagnetic 
waves by hA,hs and Ea,Eb respectively. The match- 
conditions for resonant four wave interaction are 
k E A + k E B = K A + K B > corresponding to energy- 
momentum conservation. In the center of mass system, 
all frequencies are equal and denoted by to. Thus, in this 



system we have k/^ 



-kg A , fulfilling 



the matching conditions, and corresponding to pairwise 
counterpropagating waves ^3 ■ We consider terms up to 
cubic order in the amplitudes, and note that the corre- 
sponding effects can be grouped into four categories: 

i) Nonresonant terms. Most nonlinear combinations of 
wave amplitudes are nonresonant, and their effects van- 
ishes after averaging over several wavelengths and time 
periods. 

ii) Frequency shift terms. Various combination lead to 
(small) frequency shifts of the involved waves, described 
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by interactions of the form \3Ea = ChAh* A EA + . . . . 
However, these shifts are typically not significant, and the 
corresponding influence on the conversion process will be 
neglected. 

Hi) Background curvature effects. The energy momen- 
tum pseudo-tensor of each wave generates a slowly vary- 
ing background curvature. The path of each wave will 
be described by (cubic) products of (second order) back- 
ground curvature terms and linear terms. This curva- 
ture effect on photon-graviton conversion will be small 
provided the size of all the pulses are limited in all di- 
mensions by L, fulfilling L <C kcW/oj, where W is the 
energy density of the pulses, and k is Einstein's gravita- 
tional constant. 

iv) Resonant energy conversion terms. Cubic terms in- 
volving all wave amplitudes in accordance to the match- 
ing condition correspond to pair production of photons 
by two gravitons, or the converse process. 

We will only include the energy conversion terms from 
point iv above, and the final interaction equations will 
then be of the form 

□ $1 = M$ 3 $4$2 ( 3 ) 

and permutations thereof allowed by the wave match- 
ing. Here <f>i t 2 and $3,4 represent mutually counterprop- 
agating wave fields, (i.e. Ea,b and h,A,B respectively), 
□ = d 2 — d 2 — dy — 9|, and the coefficients M depend on 
the scattering angle, the polarization of the waves and 
on the frequency, u). 

We introduce perturbations of the Minkowski met- 
ric, <?p„ = ?7p„ + h^. The perturbed energy momen- 
tum tensor T a b should then satisfy Einstein's equations 
Gat = K-T a b, where G a b is the perturbation to the Ein- 
stein tensor due to the metric perturbation. 

In order to determine the nonlinear gravitational re- 
sponse to the gravitational wave perturbations, we first 
consider the vacuum case. The perturbed metric is 

<V = w + h$ + h$ ( 4 ) 

where the hj}} are wave perturbations to first order in 

amplitude and the h^v are the second order response to 
these perturbations. For the linear wave perturbations 
we choose the transverse and traceless gauge Jl^J, such 
that the only nonzero components are hrj^ = —^22 = 
h + (t, z) and — /4i = h x (t,z). To second order, the 
vacuum Einstein equations takes the form R^,' + R~V = 
0, where R.Q is linear in /i+. x and R^? is quadratic in 
/i+,x and linear in /i/,„ . From the second order vacuum 

Einstein equations, we see that we only need the following 

(2) 

components of h)iv, corresponding to a choice of gauge: 

(2) (2) (2) (2) (2) rm. 

/ig , h 03 , h\i, h; 22 and h 33 . For the metric (gj we choose 
the tetrad basis e a vac ^ = e^^d^ — (ejj + ea^)d^ 



where 

e^ = (l + ±h®)d t -hh®dz (5) 
e[ vac) = (1- \h + - ±h<$ + \h\ + \h\)d x - \h x d v (6) 

e ( 2 vaC) = -\h x d x + (1 + \h + - |/4? + \h\ + \h 2 x )d y (7) 
et c) = h^d t + (l-lh^)d z . (8) 

In this basis the only nonzero components of R^ b 
are R$ = -R { 22 ] = 1/2(8? - d 2 z )h+ and = 
l/2(d 2 — d 2 )h x . The nonzero second order compo- 
nents are Rj^ , Rj^ , i£ x ? , R^2 and Rj^ where R^ = 
R-ab i^^ab >d 2 h+-. d 2 h 2 x ) . The second order nonlinear re- 

(2) 

sponsc terms h^J arc related to the linear wave perturba- 

(2) . . (2) (2) 

tions by R ab = 0, which gives h\i — h 22 . Furthermore, 

(2) (2) 

adding Rq to R 33 , using the approximation dt ~ ±d z 
on /i+, x , we get 

= h+d 2 h+ + h x d 2 h x . (9) 

For two (counterpropagating) gravitational waves Ka 
and h B , we have h+ = Ha+ + hs+ and h x = 
IiAx + h Bx , where the Ka and h B are given by 
Ha = tiA{t, z) exp [i(khZ — lot)] + c.c. and h B = 
h B (t, z) exp [i(— khZ — tot)] + c.c. From the form of the 
evolution equations J3J, we see that the only terms con- 
tributing to the four wave interaction are terms propor- 
tional to hAhs (and h* A h* B ), i.e. terms oc exp [— 2iujt] (and 

(2) 

exp [2iuit]). By examining R ab we see that the only com- 
ponents containing terms with this purely temporal de- 
pendence are , h^ 22 and /133 . Inserting the ansatz 
for h A and h B into Eq. ©, we get fifi — (h\ + h 2 x )/A, 

and the equation Rqq = then gives h 3 2 3 = —h$. These 
second order quantities will enter the evolution equations 
||TJ and J2J) , and contribute to the interaction coefficients 
in Q. 

The directly resonant terms described above are not 
the only contribution to the coupling coefficients. Let 
E = E A + E B = (Eai + E B i)ei + (E A 2 + E B2 )e 2 + 
(Ea3 + E B3 )e 3 and similarly for B. Furthermore, 
let Eai — E A i(t, x, z) exp [i(kEx% + ^EzZ — u>t)] + c.c, 
Ebi = E B i(t, x, z) exp [i(— k Bx x — kftzZ — u>t)] + c.c. etc. 
Inserting this into the nonlinear wave equations for the 
electromagnetic fields, the largest terms on the right hand 
of Eqs. JU and (j2J) are terms of the form Eh,Bh etc., 
which induces oscillations at frequencies and wavenum- 
bers that are nonresonant (nr), i.e. oj nr ^ k nr . These nr 
fields may combine with fields of appropriate frequency 
and wavenumber in Q and the corresponding evolution 
equations for Ha, h Bl making the resulting terms reso- 
nant with the original wave perturbation, thus contribut- 
ing to the coupling coefficients in (J3J. 

In order to take these contributions into account we 
need to modify our ansatz for the electromagnetic field 
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components, to include the nr fields. For this pur- 
pose, we use the lowest order approximation 9tB = 
— V x E, where V = (d x ,d y ,d z ), in all second order 
terms in Eq. Q). We thus make the ansatz E[ tot ^ — 

? (nr) 



separately. We use the gauge h^l^ = 0, and we choose a 
general tetrad (i.e. a basis where all are included), 

such that the extra terms corresponding to the nr fields 
basis may be added to the tetrad We only 



Eai + Ebi + E\ etc. for the total electromagnetic need to include terms up to 0(}v nl , ') in the Einstein 
field. Inserting this into Eqs. including terms 

that will contribute to the four wave coupling, we see 
that the terms linear in Iia,b takes the form DE^ nr ^ — 
ciE\hA J r c2E* A hB+CzE i ^hA+CAE* B hB J rc.c., where the Ci 
are coefficients determined by the field geometry. Thus, 
we find that the nr fields will contain oscillating factors 
exp [i(kfi - kh)} + exp [i(kg + kh)\ + c.c. Applying this 
together with the previous expressions for the original 
perturbations, noting that for the given field geometry we 
may use k Bz — ^ cos 8, k Bx = w sin 8 and — to, one ob- 
tains explicit expressions for the nr electric and magnetic 
field components. These expressions are greatly simpli- 
fied when expressing the nr field components in terms of 
the linear polarization states for the original waves. We 
thus define the two independent linear electromagnetic 
polarization states according to: Eai — — cos 9Ea+ , 
Ea2 = E A x, E A 3 = sin8E A +, and E B i = cos 8E B +, 
Eb2 = Ebxi Eb3 = — siti8Eb+- The parts of the nr 
field components contributing to the four wave interac- 
tion take the form 



tensor, since higher order terms in h^J^ will be of at 
least fourth order in the electromagnetic amplitude. The 
nonzero components of the Einstein tensor to this order 
arc the purely spatial components of G^" . By letting 

^aB ~ > ex P [ — 2iu;£] + c.c. one obtains explicit ex- 
pressions for the nonresonant gravitational fields from 
Einstein's equations. 

We now have the ingredients needed to calculate the 
evolution equations for the amplitudes of the waves and 
explicitly give the coefficients in J2J. We expand our 
metric ansatz (0J, in order to include the nonresonant 

GW fields, as = rj^ 



h TT 

IL uu 



l {2) 



+ h^v \ and in- 



clude the driven nonresonant part in our ansatz for the 
electromagnetic fields. The final tetrad basis is given by 



oJ u(nr) 



fi(vac) fi(nr) 



)dfj,, where 







E 



(nr) 



- cos 8 (— cos 8E A +h* A+ + cos 8E* B+ h B + 
E Ax h* Ax +E Bx h B x)e i( - kE -^> T +c.c. (10) 
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(11) 



and where we put = fij^ = 



m e„ Using 
this basis, and inserting explicit expressions for all fields 
into Eq. one obtains equations, involving several time 
scales, for the field components. Taking an average, over 
several wavelengths and periods, of the equations, terms 
resonant with the original perturbation are obtained |l4| . 
Expressing the equations using the linear polarization 
states, we finally obtain 

UE + = M+HiE* +±M x H n E* x , (12) 

Ax Bx B+ 

UE n+ = M+HjE* + T M x HjiE* x , (13) 

Bx Ax A+ 

from Eq. QJ, where Hi = hA+hs+ — fiAxfiBx, Hii = 
h A +h Bx +h Ax h B+ , M + = Lu 2 (l + a 2 )/2, M x = uo 2 a and 
a = cos 8. 

The evolution equations for the gravitational wave 
amplitudes are G\\ — G22 — k{Tu — T22) for h + and 
G12 + G21 = k{Ti2 + T21) for h x . Inserting explicit ex- 
pressions for all fields, we obtain the following equations 

U\h + = ±m + E I h* . +m x E n h* x (14) 

Ax Bx B+ 



plus an analogous term with spatial dependence 
exp [i(k E +k h )]. 

Terms proportional to EaEb and E A E* B oscillates as 
exp (— 2iwt) and exp (2iurf), respectively. These terms in- 
duces nonresonant GW fields h^ nr ', with the same tem- 
poral variation, through Einstein's equations. Nonlin- 
earities in the Einstein tensor of the form h^ nr ^h A B are 
then resonant with the original perturbations, and will 
contribute to the coefficients in the evolution equations 
for the GW amplitudes. 

Calculating the part of the energy momentum ten- 
sor varying as exp (±2iu;t), we make the ansatz Eai — 
EAi(t, x, z) exp [i(kE X x + k Bz z — ujt)] + c.c. etc. We have 
not included the nonresonant fields calculated above in 
the ansatz, since these fields do not appear to second 
order. The nonzero components of the energy momen- 
tum tensor varying as exp (±2iwi) are the purely spatial 
components of T a b. 

Since the energy-momentum tensor induces grav- 
itational fields with the same temporal variation, 
exp(±2iwi), we assume a general metric perturbation 
with this time dependence. We can separate out terms 
in Einstein's equations describing the GWs and the non- 
linear response to these, and terms describing the metric 
response to the energy momentum tensor of the EMWs, 
and therefore we can calculate the nr gravitational fields 



Dh 



±m+E I K k + +m x E H h* x (15) 

Ax A+ 



where £7 = E A +E B + + E A xE B x, E u = E A +E B x ~ 
EaxEb+, m+ = k{1 + a 2 ) and m x = 2na. The coupled 
system of equations <|12D (I15|l describes the evolution of 
the slowly varying wave amplitudes. 
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Introducing the energy densities of the waves Saw ~ 
co 2 /2K(\h+\ 2 + \h x \ 2 ) and Sem = \E+\ 2 + \E x \ 2 for the 
gravitational and electromagnetic wave respectively, we 
note that the total wave energy, Etot-, is conserved in the 
interaction described by Eqs. I|12f> -H15 [) ■ 

In order to illustrate the physics of |Q — we as- 
sume long pulses, i.e. □ — * — 2iw9j. For simplicity, we 
consider the case where all waves have +-polarization, 
and assume that the initial amplitudes and relative 
phases of the wave envelopes are fixed in such a way 



that Ea- 



E 



B + 



E and h 



A+ 



1_B + 



h where 



E = E(t)exp[iip E (t)] and h = h(t)exp[i(p h (t)]. We 
insert this into Eqs. (|12|) and (|14|l . and reexpress the 
equations as evolution equations for the energy den- 
sities Sem and Sow, an d the relative phases of the 
waves. Then, for the special choices of initial values 
(iph — ifE)\t=o — ± 7r /4, one obtains very simple analyti- 
cal solutions for the ratio of electromagnetic and gravita- 
tional wave energy densities. The results in this case are 
£gw/£em = (£aw/£EM)(0)e~ T in the case when (<p h - 
<Pe)(0) = -7r/4 and £em/£gw = {£em /£aw)(0)e~ T for 
{ip h - ip E )(0) = 7r/4, where r = (1 + a 2 )K£ tot t/2uj. 

Thus we see that the characteristic time-scale for 
graviton-photon conversion is rather slow. In particu- 
lar, as far as the interaction between four long pulses 
is concerned, there is only time for large energy conver- 
sion if the waves also contains enough energy to signif- 
icantly modify the background curvature (c.f. point Hi 
above). However, we note that the GW-EMW conver- 
sion takes place in a space-time that is locally flat, and 
a given pulse (EMW or GW) can interact with other 
pulses (or even individual quantas) consecutively, leading 
to significant conversion after multiple processes. Thus, 
we here emphasize that the coefficients in Eqs. (11211 l|15|) 
can be used also to describe incoherent interaction involv- 
ing multiple processes. Writing the equations in terms 
of the vector potential, A, rather than the electric field 
and making a normalization h — > h/^/2n, we note that 
M + ,m+ — ► M+ = klu 2 (1 + a 2 ) and M x ,to x — > M x = 
2nuj 2 a in lfT2 )l -l(T5 |l . The coefficients M+ and M x now 
correspond to the scattering amplitude matrix elements, 
from which the cross-section can be calculated straight- 
forwardly. Following Ref. [l^, we find that the differen- 
tial cross-section for graviton-photon conversion averaged 
over all polarization states is 



da 



\M+\ 2 + \M, 
128w 2 (2tt) 2 



(l + « 2 ) 5 



(16) 



where Lp is the Planck length. 

The process of graviton-photon conversion may lead 
to thermalization between gravitons and photons in the 
early universe. For thermalization to occur before ex- 
pansion slows down the conversion rate too much, the 
collision frequency must fulfill v = an > H, where H 
is the Hubble parameter and n is the number density. 



Considering a photon gas at temperature T in the ra- 
diation dominated era, using H ~ 1/t, the above con- 
dition can be written T/T P > W c ^ 6 where T P is the 
Planck temperature and Af e s is the effective number of 
degrees of freedom. Thus thermalization, in case this is 
still a meaningful concept when T ~ Tp, must be de- 
scribed within a quantum theory of gravity. The classi- 
cal photon-graviton conversion process nevertheless puts 
some limits on the background levels of gravitational ra- 
diation, since it makes a relic GW spectrum much above 
the thermal level unlikely. 

We note, however, that the time-scale for coherent in- 
teraction is shorter than that for incoherent single parti- 
cle interactions by a factor of the order k 2 L 2 P: and thus 
physical scenarios with significant GW-EMW conversion 
in astrophysics is still possible within a classical frame- 
work. Due to the large energy densities required, how- 
ever, such processes are likely to take place in signifi- 
cantly curved backgrounds, for which a description lies 
outside the scope of this letter. 

A comparison of the cross-section (lltifl to that obtained 
by using Feynman diagrams would be of much interest, 
but as far as we know, the later result is yet to be de- 
duced. We note, however, that an interaction lagrangian 
of the form C cx WT^, as used in e.g. Ref. [llj in lin- 
ear quantum gravity, will not be sufficient here, since the 
effect considered in this letter is quadratic in the gravi- 
tational wave amplitudes. 

Andreas Kallberg, who was supported by the Swedish 
National Graduate School of Space Technology, was trag- 
ically killed in an accident shortly after this work was 
completed. 
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